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In this short note we translate the best available observational bounds on the CMB bispectrum 
amplitudes into constraints on a specific scale-invariant New Physics Hypersurface (NPH) model of 
vacuum state modifications, as first proposed by Danielsson, in general models of single-field infla- 
tion. As compared to the power spectrum the bispectrum constraints are less ambiguous and provide 
an interesting upper bound on the cut-off scale in general models of single-field inflation with a small 
speed of sound. This upper bound is incompatible with the power spectrum constraint for most of 
the parameter domain, leaving very little room for minimal cut-off vacuum state modifications in 
general single-field models with a small speed of sound. 

INTRODUCTION 

Bounds on higher-order non-Gaussian contributions to the observed and approximately Gaussian temperature 
anisotropy map of the Cosmic Microwave Background are rapidly entering a phase where they can be used to constrain 
and sometimes even rule out certain (exotic) models of inflation. For instance, bounds on the amplitude of the CMB 
three-point function, or bispectrum, have already narrowed down the parameter range of string inspired Dirac-Born 
Infcld models of inflation. 

In a separate development it was recently shown that modifications to the inflationary vacuum state can also result 
in significantly enhanced non-Gaussian signals [IH1]- In particular the bispectrum can reveal enhanced features that 
essentially arise due to interactions between excited quanta in the modified vacuum. The power spectrum is also 
sensitive to deviations from the Bunch-Davies vacuum, but there the corrections are suppressed by the Bogoliubov 
parameter, whose magnitude is small and typically related to the inflationary Hubble parameter divided by the scale 
of new physics, which is already constrained to be smaller than 10 -2 [5]. For the bispectrum one instead encounters 
an additional enhancement factor, independent of the Bogoliubov parameter, that is in fact proportional to the scale 
where effective field theory is supposed to break down. The appearance of this enhancement factor can be traced 
back to the New Physics Hypersurface (NPH) used to define the modified vacuum state. The higher the cut-off 
scale representing new physics, the more time the fields have before they cross the horizon and freeze, and the more 
time the interactions effectively have to generate a significant non-Gaussian contribution. This also implies that the 
(decoupling) limit of taking the cut-off scale to infinity in a modified vacuum is in general ill-defined, producing an 
infinite bispectrum, whereas the power spectrum yields the standard result in this decoupling limit (assuming the 
Bogoliubov parameter is suppressed by one over the cut-off scale). 

In this work the most recent bispectrum constraints will be applied to a particular and arguably the simplest 
model of vacuum state corrections, which was first considered in [6], but constructed and applied more generally by 
Danielsson [7] . We will derive to what extent this vacuum state model can be ruled out by the bispectrum data alone 
in two particular scenarios; a slow-roll single field inflationary model with a higher derivative operator representing 
the interactions and a non-canonical single field model with a small speed of sound. The bispectrum prediction for 
a general (NPH) vacuum state modification in these single-field models was computed in [3l H] . As an important 
and hopefully illuminating illustration we will slightly extend and improve the analysis, and apply it to the specific 
vacuum state proposal of Danielsson. 

This paper is organized as follows. We will start with a quick review and small generalization of the Danielsson 
vacuum state proposal, followed by a summary of the results for the bispectrum in the context of a vacuum state 
modification in general single-field models of inflation. We then explain how these results can be turned into an upper 
bound on the scale of new physics and to what extent this vacuum state proposal is ruled out or not. We end with 
some conclusions and prospects for future improvement of the bounds. 
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A MINIMAL MODEL OF VACUUM STATE MODIFICATIONS 



Let us begin with a short review of Danielsson's original proposal for a modified vacuum state, under the natural 
assumption that a high energy cut-off scale L c exists beyond which effective field theory breaks down [7]- Solving for 
the mode functions in an inflating (pure de Sitter) background one writes down solutions for the field and conjugate 
momentum operators in terms of creation and annihilation operators in the standard way. Now identify a (conformal) 
time rjo serving as the initial time where one would like to define the vacuum state. This initial time 770 is related to 
the cut-off A c by demanding that the physical momentum at 770 equals the cut-off scale, i.e. \krjo\ = -jj-. This ensures 
that the description is always within the effective field theory regime. Note that since we are imposing a cut-off on 
the physical momentum the initial time 770 is necessarily a function of the comoving momentum k. With respect to 
this initial time 770 the creation operators (and subsequently the annihilation operators) can then be expressed as 

a k (v) = Ukiv) a k (vo) + Vk(v)aL k {vo) , (1) 

describing nothing else but the mixing of creation and annihilation operators as time progresses. To ensure the time- 
independent commutation relations for the creation- and annihilation operators requires that |7ifc| 2 — \v k \ 2 — 1. Let 
us now define a natural candidate for a vacuum state at some initial conformal time r)a 

^(770)10,770) =0. (2) 

Obviously this choice requires that v k (ri) vanishes at 77 = 770 and as a direct consequence this candidate vacuum 
corresponds to a minimal uncertainty state at 77 = 770 [7 . This choice of vacuum can be understood as selecting the 
'local' empty state at the time the physical momentum equals the high-energy cut-off A c . This time will be different 
for different comoving momenta and defines a so-called 'New Physics Hypersurface'. It can therefore be viewed as the 
vacuum state that is closest to the Bunch-Davies state in the presence of a high-energy cut-off and can in that sense 
be regarded as a minimal modification. 

To understand the relation of this choice of vacuum to the standard Bunch-Davies state, consider the field operator 
which should be proportional to the sum of the creation and annihilation operators. In terms of mode-function 
solutions to the equations of motion f k (77) this implies that 

fk(v)=N k (u k (r,) + vt(v)) , (3) 

where the overall (real) normalization N k is fixed by the Klein-Gordon normalization condition. Similarly, the expres- 
sion for the canonical momentum operator defines a function g k (r)) 7 which is proportional to the difference between 
1^(77) and v k (r))- 

g k {v) = N k (ufcfa) - v%(Tj)) . (4) 
Using these relations one can derive the following expression for the function v k (77) in terms of f k and g k 

w£fa) = \ K" 1 Mv) ~ N- 1 g k ( V )] . (5) 

Defining a vacuum state a k (r)o) |0,?7o) = is now explicitly seen to be equivalent to picking a mode-function solution 
f k (and its canonically conjugate function g k ) such that v k (rf) vanishes at some 77 = 770. Bunch-Davies mode-functions 
have the property that 7^(77) only vanishes in the limit 7/ — > —00, i.e. the Bunch-Davies state in this class of vacua 
corresponds to the minimal uncertainty vacuum state in the infinite past. 

Now let us instead consider the general solution, which can be constructed from the Bunch-Davies mode-function 
and its complex conjugate 

fk(v) = A k f k BD (v) + B k f k BD *(v), (6) 

where A k and B k are complex coefficients satisfying \A k \ 2 — \B k \ 2 — 1. Similarly, the momentum mode-function g k (rj) 
can be expressed in terms of the Bunch-Davies momentum mode-functions 

g k (r ) )=A k g* D (r 1 )-B k gZ D *( V ). (7) 

To find the natural candidate vacua at some fixed time 770 one demands that v k (rjo) = 0. Using ^ this gives the 
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following expression for the Bogoliubov rotation parameter fa = ^ in terms of the Bunch-Davies field and momentum 
mode-functions 



In terms of the (Bunch-Davies) functions Uk BD (j]) and Vk BD (rj) that appear in the expression for the annihilation 
and creation operators ([I]) the result simply reads 



u l Wo) 



A couple of comments are in order. Since the (late time) power spectrum of fluctuations is proportional to \fk\ 2 , 
the Bogoliubov rotation parameter fa denotes the leading order correction to the power spectrum due to a vacuum 
state modification. Scale invariant vacuum states of course require that fa is independent of the comoving momentum 
k, which will be guaranteed by the relation krj = — % , i.e. 770 is fc-dependent. Also note that this minimal cut-off 
vacuum construction has so far been completely general, not depending on any details of the inflationary Lagrangian. 
Whatever inflationary model one is interested in, to construct these states and determine the Bogoliubov parameter 
fa one simply plugs in the relevant (mode-) functions in the inflationary case of interest. 

The best known (and originally discussed) example is that of the massless scalar field, that is associated to standard 
slow-roll inflation. The Bunch-Davies field and momentum mode-function in that case read 

A) = ^'H) ' ^ D (v) = ] fle~^. (10) 

"k 



Using the expression Q or ^ one then derives the following expression for the Bogoliubov parameter ft, 



fa = ^ ± —^ 2lkll ° • (11) 

Since \kr)o\ = ^> 1 this Bogoliubov parameter is indeed scale invariant and approximates to fa ~ -M-e^i* - ^. 

As we have emphasized, the construction of these minimal cut-off states is completely general and can just as well 
be applied in the context of small speed of sound models of inflation, which includes DBI inflation as a special class. 
One might suggest that in non-canonical models of inflation the introduction of a cut-off scale is not required, since an 
infinite number of higher order corrections have been resummed leading to the non-canonical kinetic term. However, 
these corrections are only a subset of all the possible higher-dimensional operators that are expected to contribute as 
the string- (or cut-off) scale is approached, motivating the general introduction of a cut-off scale 1 . So the main, and 
for our purposes only, difference between small speed of sound and slow-roll models of inflation is the introduction of 
a reduced (and assumed constant) speed of sound c s < 1 of the inflaton fluctuations. For the Bunch-Davies functions 
u BD and v BD one can check that the only effect is the appearance of a single factor of the speed of sound c s in the 
krjo terms as compared to the standard slow-roll case. The end-result for the Bogoliubov parameter therefore reads 

fa = — % - e - 2ikc ^° , (12) 

2fcc s ?7o + i 

which translates into the scale-invariant approximate expression fa 2 B A e'(s g h~ ). One observes that the ratio 
H / A c governing the magnitude of the corrections in standard slow-roll is modified to H/c s A c in these models, which 
could be orders of magnitude bigger. Note that the length scale 1/H* = c s /H < 1/H corresponds to the sound horizon 
in small speed of sound models, corresponding to the length scale at which the behavior of the mode-functions turns 
non-adiabatic. The appearance of the factor kc s r) is a general feature in these models, which was for instance also 
observed for the enhancement factors in the three-point function that appear whenever one introduces an arbitrary 
modified initial state at some high-energy cut-off. We conclude that the same is true for the Bogoliubov parameter in 



1 More specifically, in [8] it was shown that in DBI models of inflation under certain conditions the standard procedure for defining the 
Bunch-Davies vacuum can fail. 
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cut-off modified initial state proposals: the magnitude (and phase) depends on the combination kc s rjo which results in 
larger absolute values of the Bogoliubov parameter for the same value of the cut-off and inflationary Hubble parameter, 
as compared to standard slow-roll models. 

This ends our short summary of minimal cut-off initial states. The results for the three-point function, as reported 
in [31 3] , were model-independent in the sense that they only relied on the presence of a physical momentum cut-off. 
The Bogoliubov parameter was left as a free parameter and general bounds were derived on its magnitude for different 
inflationary Lagrangians. In this work we would like to specifically constrain the minimal cut-off modified vacuum 
state models for which we will need the expressions for the Bogoliubov parameters derived above. 



As already referred to, in our previous work we calculated the results for the three-point function under the 
assumption of an arbitrary scale-invariant initial state modification in the NPH framework. The presence of a 
physical high-energy momentum cut-off was assumed, which identified an initial time for each comoving momentum 
mode where the initial state was defined by introducing a k- independent, but undetermined, Bogoliubov parameter. 
The presence of such a cut-off alone already has major consequences for the three-point function, also known as the 
bispectrum in Fourier space. Whenever the Bogoliubov parameter is non-zero, particles are injected at the initial 
(cut-off) time, which allows potential (self-) interactions to generate a large bispectrum amplitude at the time of 
horizon crossing. In the calculations these effects show up as large enhancement factors that depend on powers of 
A c /H 3> 1. Depending on the details of the inflationary model under consideration the power of the enhancement 
factor could be as large as three, resulting in a huge bispectrum amplitude. The only reason why these large non- 
Gaussian signals have not yet been detected or ruled out is their specific (oscillatory) shape, which is orthogonal to 
any of the non-Gaussian templates used in analyzing the CMB data so far. Proposals for better adapted templates 
and improved methods to detect or constrain these oscillatory non-Gaussian signals have recently been put forward 
[9l IIP), but have not yet been applied to the available data, forcing us to concentrate on the standard non-Gaussian 
shapes that have been constrained. Projecting the large bispectrum onto any one of the available templates still allows 
one to derive reasonable constraints in general. In fact, in some cases the constraints are already stronger than those 
available from the power spectrum. Before discussing the results of the projections to the observational templates 
it should be pointed out that in principle the results only apply to the three-dimensional bispectrum. However, in 
general the changes resulting from the reduction to the two-dimensional sphere are minimal and can (partially) be 
taken into account by introducing a weight function in the three-dimensional analysis 

For models with a small speed of sound (c s <C 1) the projections to the local and orthogonal templates [12], including 
the phase, were computed analytically [5j . To derive these results the assumption was made that the modifications to 
the BD state were strictly oriented in the k\ direction, introducing a constant kirjo parameter, breaking the symmetry 
between the different momenta in the triangle. Properly maintaining the symmetry would instead require introducing 
three constant parameters kitjo, where fcj is the direction in which one has perturbed the BD vacuum state via a 
Bogolyubov transformation. Consequently, one can factor out the constant fc^o for each perturbed state ^ . As 
expected, the results of this more complete analysis that preserves the symmetry between triangle momenta differ 
only slightly. Moreover, we have ameliorated our (numerical) integration methods and used the improved definition 
of the inner product between shapes as proposed by which introduces a weight function to simulate projection 
effects onto the 2D CMB sky. This increases the correlation between smooth and oscillating shapes as oscillations are 
slightly damped by the weight function. 

General single field models are defined by a single (Lagrangian) function P(X, 0), where X = d^4>d u 4>. From this 
function P(X, <j)) expressions can be derived for the slow-roll parameters, the speed of sound and two other variables 
E and A, whose ratio ^ will appear in the three-point function [TJEI]. Introducing a modified NPH vacuum state 
the final enhancement factor in the projections turns out to be at best quadratic in p = kc s r]o = A c c s /H 1. The 
results, for an undetermined Bogoliubov parameter j3 = \(3\e lS , in the limit that \ ^> ^ and p ^> 1 is 



GENERAL BISPECTRUM PREDICTIONS 



A j*local 

^JNL 



a/; 



4 x 10 _ 
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cos (S), 



cos(<5) . 



(13) 



(14) 



The additional terms that have been neglected are subleading in p and are out of phase in the sense that the phase 
parameter 5 is shifted with tt/2. Also note that these results do not reduce to the standard slow-roll results in 
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the c s = 1 limit. Subdominant terms that have been neglected in the c s -C 1 limit will turn in to the dominant 
contributions to the bispectrum in the c s = 1 limit, explaining why the c s = 1 limit of the above result does not 
reproduce the slow-roll bispectrum. It is also important to stress that the DBI models are a special class of small 
sound speed models for which 2A/S = 1 — l/c 2 s and as a result the leading contributions to the bispectrum cancel 
exactly [13] • The remaining contributions to the projections of the DBI modified NPH initial state bispectrum are 
only linearly enhanced in the small c s limit and the details of the projections in the DBI case can only be evaluated 
semi-analytically [3]. This leads to DBI constraints that are far less interesting as compared to general small sound 
speed models. We will therefore concentrate on the general models and not discuss the DBI results in detail. 

In the case of standard canonical single-field slow-roll inflation with a specific dimension 8 higher-derivative term 
[14] the enhancement factor appearing in the bispectrum amplitude was found to be quadratic in A c /H [3 . However, 
in the projection to the available observational templates f° L = A(0, T) f^ L , where A(T, O) is the projection factor, 
one looses one power of the enhancement, leaving a linear enhancement in A c /H as far as the projections are concerned. 
In addition the amplitude of the enhancement is further reduced by slow-roll (via e) and the dimensionful coupling 
of the higher-derivative term. In these results the effects of a (cut-off dependent) phase in the Bogoliubov parameter 
was neglected and for our purposes here we would like to remedy this situation. To incorporate the phase we had to 
rely on numerical methods to determine the projection factors. For the projection onto local non-Gaussianities we 
obtain (in the limit of large p = A c /H) 

A/^T 1 * |pAe|/?| (^^co S (5), (15) 

which, not surprisingly, is very similar to equation (6.11) in [3J, except for the appearance of a phase dependence. 

It is important to realize that the inclusion of the perturbative effect of the dimension 8 higher-derivitive operator 
is only consistent when </>/A 2 < 1 [TJ], which translates into e > A^/M^H 2 . The power spectrum observations tell 

us that e = which implies that H/A c < 10~ 2 to ensure a consistent higher-derivative expansion. As a 

' pi 

consequence the non-Gaussianities due to the presence of the dimension 8 higher-derivative term reach a maximal 
magnitude for H/A c ~ 10 -2 , yielding an equilateral amplitude of at most f^ L ~ 0(1). On the other hand, modifying 
the initial state leads to an additional enhancement factor that is proportional to the ratio between the UV cutoff A c 
and the scale of inflation H . The required consistency of the higher-derivative expansion, being inversely proportional 
to the cut-off scale A c , reduces this enhancement and as a consequence weakens the constraints on /3. In fact, as 
we will see in the next section, for the minimal cut-off vacuum state proposal assumed in this paper the constraint 
becomes practically meaningless. This should not come as much of a surprise, since in order to boost non-Gaussianities 
from higher-derivative operators the cut-off scale needs to be relatively close to the inflationary scale, whereas the 
non-Gaussian amplitude as a consequence of initial state modifications grows as the the hierarchy between the cut-off 
and the Hubble scale increases. For completeness, let us remark that the correlation between the orthogonal template 
and the standard single-field higher-derivative bispectrum, being very similar, yields comparable constraints. 

After this very short review and slight improvement of what is known about bispectrum projections due to a generic 
NPH initial state, let us now derive the best available constraints on the parameters of minimal cut-off initial states. 



BISPECTRUM CONSTRAINTS ON MINIMAL CUT-OFF VACUUM STATES 

Having summarized the results for minimal cut-off initial states and the general projected non-Gaussian amplitudes 
due to arbitrary cut-off (or NPH) modified vacua, we would now like to apply the most recent observational bounds 
to derive constraints on the parameters of the minimal cut-off initial state proposal. 

To proceed we will make use of the following recent (2 a) constraints on an d /nl 031 US] 

-10 </;° L cal < 74, 
-410 </S r L *< 6 

The observational constraints on the equilateral amplitude /S^ 1 can be ignored, because the projection to the 
equilateral template turns out to be orders of magnitude smaller. One should keep in mind that the power spectrum 
constraint on |/3| is roughly 10 -2 , which is based on the lack of evidence for oscillatory behavior with a larger amplitude 
in the power spectrum [5]. 

A notable feature of the results for the projections is that they oscillate as a function of the phase 5 of the Bogoliubov 
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parameter ( 13 14 ). Different values for the phase can therefore result in constraints that deviate considerably. Indeed, 
for special values of the phase the projection vanishes and the constraints disappear altogether 2 . In minimal cut-off 
models the phase of the Bogoliubov parameter is actually a function of the cut-off scale A c (11-121 and as a consequence 



for special values of A c the projected amplitudes would be significantly reduced, as can be seen in figure ([!]). Since 
the precise value of the cut-off A c is unknown, the constraints will be derived under the reasonable assumption that 
a typical value of A c is expected to be closer to an (absolute value) maximum in the projection than to the special 
point where the projection vanishes. To implement this we will derive the constraints based on an expression that is 
related to the absolute value average over a single oscillatory domain in the parameter A c . 
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Figure 1: Contribution of initial state modifications to orthogonal non-Gaussianities for various values of p/c 2 .The red upper 
and lower parts in the diagram correspond to the observationally excluded regions. 



Let us first analyze general single field models with a small speed of sound. The minimal cut-off vacuum state 
proposal predicts a Bogoliubov parameter with an absolute value approximately equal to \/3k\ ~ 2< /a an< ^ a P nase 
5 — i(§7T — )■ In the limit p = A c c s /H > 1 we derive, to leading order in p, 



A -/"local 
^JNL 



2 x 1CT 3 4 sin(», 



A/° r * ~ 3 x 1CT 2 ^ sin(2p), 
ci 



(16) 
(17) 



The absolute value of the Bogoliubov coefficient has removed one enhancement factor, resulting in a linear p depen- 
dence of the prefactor in front of the oscillatory sin(2p) term. As an example, the maximum contribution to the 
orthogonal non-Gaussian template from initial state modifications has been plotted in figure [l] As already alluded 
to, the oscillations in the projection factor force us to work with the average over a single oscillatory domain in order 
to derive a constraint on the ratio A c /H. Simply replacing the sin(2p) with the average 2/7T the derived bounds 
would only apply to about 60 percent of the p domain. If instead demanding that at least 90 percent of the domain 
is included within the constraints, we need to replace sin(2p) with ~ 0.16. 

The ambiguous sign due to the oscillatory behavior also prompts us to consider the bound on the (largest) absolute 
value of the non-Gaussian amplitude, which for the local amplitude is |/]^£ al | < 74 and for the orthogonal type 
\Inl\ — 410- Using these largest absolute values for the local and orthogonal amplitude and insisting on a 90 percent 



2 To be more precise, for those values one can no longer neglect the subleading terms in p, which are in effect tt/2 out of phase and 
therefore maximize exactly when the leading terms minimize. Nevertheless, the fact that they are subleading in p implies that the 
constraints for those values of p become practically meaningless. 
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coverage, an upper bound on the parameter y£ can De derived from local and orthogonal type non-Gaussianities 

local: — <2.5xl0 5 c s , orthogonal: — <8.5xl0 4 c s (18) 
H H 

Clearly the strongest result is obtained from the orthogonal bound, becoming stronger as c s is smaller. In fact, 
the speed of sound c s cannot be much smaller than 10 -2 because that would produce a leading order equilateral 
non-Gaussian result in conflict with observations. Applying that minimal value of c s the upper bound on 4f roughly 
equals 

^ < 10 3 , (19) 

Again we would like to stress that this is an upper bound. The cut-off scale should be close to the inflationary Hubble 
scale to avoid too large an enhancement of the (projected) non-Gaussian signal, with the exception of special higher 
values of A c for which the phase <5 is close to an integer number times n, corresponding to 10 percent of the parameter 
domain. This constraint is weaker than initially anticipated in 4 , mainly because the phase in this case is a function 
of A c , producing oscillations in the non-Gaussian signal, forcing us to introduce a factor related to the average of 
the absolute non-Gaussian signal and the most conservative (upper) bound on the absolute value of the observed 
constraints. 

Incorporating the, admittedly crude, lower bound on A^ > iq 2 from analysis of the power spectrum, for a general 
speed of sound one arrives at 

ID 2 A 

— < — < 8.5 x 10 4 c s (20) 
c s H 

There is a small window remaining for c s ~ 10 , but for small speed of sound models with c s ~ 10 -2 one can conclude 
that minimal initial state modifications are practically ruled out, except for those special values of the phase at which 
the contribution to the orthogonal non-Gaussian template (nearly) vanishes, corresponding to about 10 percent of 
the parameter domain. 

To briefly illustrate the impact of the oscillatory behavior in this analysis, consider neglecting the oscillating 
behavior, choosing sin2p = 1 for local non-Gaussianities and sin2p = — 1 for orthogonal non-Gaussianities. Again 
this would yield a strongest upper bound from the orthogonal constraints that is roughly equal to A c /H < 2 x 10 2 c s , 
which for c s ~ 10~ 2 can essentially be ruled out from the non-Gaussian data alone and would in general be inconsistent 
with the power spectrum data that shows no evidence of an oscillatory effect with an amplitude that large. 

As already commented on, for DBI models the constraints weaken considerably because the enhancement of the 
non-Gaussian signal is reduced by a factor of p. This results in no enhancement at all in the projection to the local (or 
orthogonal) template, removing all dependence on the cut-off scale. As a a consequence the observational constraints 
do not result in (interesting) bounds on the cut-off scale A c , which parameterizes the initial state modification. 

Finally, let us briefly consider standard slow-roll inflation (c s = 1), including a dimension 8 higher derivative term 
that, together with the initial state modification, is responsible for a large p 2 enhanced non-Gaussian signal. Inserting 



and the phase into eq. ( 14 1 one obtains 



fSt * -^sin(2p) (21 i 



All enhancement in p = A c /_ff is lost and one is left with a maximum contribution to the local non-Gaussian template 
that is of (less than) order 1, which is significantly below set constraints on local type non-Gaussianities. As a 
consequence minimal initial state modifications in this particular canonical single field inflationary model are not 
constrained by the local non-Gaussian projection. In order to become sensitive to the intrinsically large non-Gaussian 
signal set by the ratio A c /H, templates better adapted to the theoretical non-Gaussian shape should be designed and 
compared to the available data. 



CONCLUSIONS 



To summarize, we have applied the currently best available constraints on local, orthogonal and equilateral shape 
non-Gaussianities to derive bounds on the cut-off scale parameterizing a (slightly generalized) minimal model of initial 



state modifications. The results strongly depend on the model of inflation under consideration. Interesting bounds 
can be derived in the context of general, non-DBI, single field small speed of sound models of inflation. In that case, 
non-Gaussian constraints from the orthogonal template alone lead to the following upper bound on the cut-off scale 

^<8.5xl0 4 c s . (22) 
H 

Combined with results from the equilateral template and the power spectrum leaves only a very small window 
(c s ~ 10 _1 , A c /H ~ 10 3 ), and minimal initial state modifications in general single field models with a small speed 
of sound can almost be ruled out. To be more precise, the above upper bound is valid for 90 percent of the cut-off 
parameter domain, excluding small isolated regions covering in total 10 percent of the parameter domain where the 
projection to the orthogonal template nearly vanishes as a consequence of the oscillating nature of the projection. 

For DBI and a canonical single field model with a dimension 8 higher derivative operator the analysis does not lead 
to interesting non-Gaussian constraints. The reason is the complete lack of sensitivity of the available templates to 
the large non-Gaussian signal of minimal initial state modifications. Consequently, more suitable templates have to be 
developed to efficiently probe (minimal) initial state modifications in these models, perhaps along the lines of [91 110). 
Work in this direction is ongoing and we hope to report on more effective non-Gausian templates and observational 
strategies in the near future. 

An important property that we have emphasized is that (minimal) initial state modifications lead to a non-Gaussian 
signal that increases as the separation between the cut-off scale and the inflationary Hubble parameter increases. As a 
consequence non-Gaussian constraints can only give rise to an upper bound on the cut-off scale parameterizing initial 
state modifications arising due to new physics. Together with lower bounds on the cut-off scale from other consider- 
ations (including the power spectrum data or generic effective field theory arguments) this carries a strong potential 
to rule out these types of initial state modifications in the future as the non-Gaussian analysis and corresponding 
constraints improve. 

At the same time the interesting feature of an upper bound on the cut-off scale poses somewhat of a theoretical 
conundrum, since it is obviously in conflict with the idea of decoupling. Perhaps this should be interpreted as signalling 
a fundamental flaw of these type of initial state modificatitions that might be responsible for inconsistencies in the 
perturbative expansion of the quantum field theory, conceivably ruling out these vacuum states. We hope to come 
back to this important issue in future work. 
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